Asymptomatic Plasmodium falciparum carriage at the end of the dry season is associated with subsequent infection and clinical malaria in Eastern Gambia.

BACKGROUND: Chronic carriage of asymptomatic low-density Plasmodium falciparum parasitaemia in the dry season may support maintenance of acquired immunity that protects against clinical malaria. However, the relationship between chronic low-density infections and subsequent risk of clinical malaria episodes remains unclear. METHODS: In a 2-years study (December 2014 to December 2016) in eastern Gambia, nine cross-sectional surveys using molecular parasite detection were performed in the dry and wet season. During the 2016 malaria transmission season, passive case detection identified episodes of clinical malaria. RESULTS: Among the 5256 samples collected, 444 (8.4%) were positive for P. falciparum. A multivariate model identified village of residence, male sex, age ≥ 5 years old, anaemia, and fever as independent factors associated with P. falciparum parasite carriage. Infections did not cluster over time within the same households or recurred among neighbouring households. Asymptomatic parasite carriage at the end of dry season was associated with a higher risk of infection (Hazard Ratio, HR = 3.0, p < 0.0001) and clinical malaria (HR = 1.561, p = 0.057) during the following transmission season. Age and village of residence were additional predictors of infection and clinical malaria during the transmission season. CONCLUSION: Chronic parasite carriage during the dry season is associated with an increased risk of malaria infection and clinical malaria. It is unclear whether this is due to environmental exposure or to other factors.

Genomic variation during culture adaptation of genetically complex Plasmodium falciparum clinical isolates.

Experimental studies on the biology of malaria parasites have mostly been based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of some Plasmodium falciparum clinical isolates in analyses focusing on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections, which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points over several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and time points. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity, although proportions of genotypes varied over time. Drug resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (including AP2-HS, EPAC and SRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identified de novo variants not detected in the bulk isolate sequences. Interestingly, several of these were nonsense mutants and frameshifts disrupting the coding sequence of EPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of genomic identity by descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within endemic populations.

Pf7: an open dataset of Plasmodium falciparum genome variation in 20,000 worldwide samples.

We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network.  It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented.  For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations.  We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent.  We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines.  Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website.

Genetic Diversity of Plasmodium falciparum and Distribution of Antimalarial Drug Resistance Mutations in Symptomatic and Asymptomatic Infections.

Malaria control relies on passive case detection, and this strategy fails detecting asymptomatic infections. In addition, infections in endemic areas harbor multiple parasite genotypes that could affect case management and malaria epidemiology. Here, we performed AmpSeq genotyping to capture polymorphisms associated with antimalarial resistance and the genetic diversity within natural Plasmodium falciparum infections. Known genetic polymorphisms associated with altered drug susceptibility were screened for the five most common marker genes, pfdhfr, pfdhps, pfmdr1, pfcrt, and pfK13, and genetic diversity was established from two known AmpSeq markers, cpmp and csp. Relative abundance of the different genotypes within mixed infections was calculated from the number of reads per genotype. Genotyping was performed on 117 samples, 63 from asymptomatic and 54 from symptomatic individuals. We identified up to 15 genotypes within an infection, and the median multiplicity of infection was higher in asymptomatic infections (median MOI = 5 in asymptomatics versus median MOI = 2 in symptomatics, P < 0.001). No genetic differentiation on parasites from asymptomatic and symptomatic individuals was found. No mutation associated with ART resistance was identified. Prevalence of the P. falciparum chloroquine resistance wild-type genotype (CVMNK) reached 80%, confirming a return to chloroquine (CQ) sensitive parasites in Cameroon. In addition, the CQ-associated resistant genotype (CVIET) was present at very low density in polyclonal infections. Persistence of low-density chloroquine resistant parasites indicates competition-survival trade-offs may contribute to maintaining genetic diversity in natura. Thus, monitoring the expansion of these low-density genotypes in different immune backgrounds will be critical to evaluate drug policy changes.

A Cohort Study on the Duration of Plasmodium falciparum Infections During the Dry Season in The Gambia.

BACKGROUND: In areas where Plasmodium falciparum malaria is seasonal, a dry season reservoir of blood-stage infection is essential for initiating transmission during the following wet season. METHODS: In The Gambia, a cohort of 42 individuals with quantitative polymerase chain reaction-positive P falciparum infections at the end of the transmission season (December) were followed monthly until the end of the dry season (May) to evaluate infection persistence. The influence of human host and parasitological factors was investigated. RESULTS: A large proportion of individuals infected at the end of the wet season had detectable infections until the end of the dry season (40.0%; 16 of 40). At the start of the dry season, the majority of these persistent infections (82%) had parasite densities >10 p/µL compared to only 5.9% of short-lived infections. Persistent infections (59%) were also more likely to be multiclonal than short-lived infections (5.9%) and were associated with individuals having higher levels of P falciparum-specific antibodies (P = .02). CONCLUSIONS: Asymptomatic persistent infections were multiclonal with higher parasite densities at the beginning of the dry season. Screening and treating asymptomatic infections during the dry season may reduce the human reservoir of malaria responsible for initiating transmission in the wet season.

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination.

Understanding Host-Pathogen-Vector Interactions with Chronic Asymptomatic Malaria Infections.

The last malaria parasite standing will display effective adaptations to selective forces. While substantial progress has been made in reducing malaria mortality, eradication will require elimination of all Plasmodium parasites, including those in asymptomatic infections. These typically chronic, low-density infections are difficult to detect, yet can persist for months. We argue that asymptomatic infection is the parasite's best asset for survival but it can be exploited if studied as a new model for host-pathogen-vector interactions. Regular sampling from cohorts of asymptomatic individuals can provide a means to investigate continuous parasite development within its natural host. State-of-the-art techniques can now be applied to such infections. This approach may reveal key molecular drivers of chronic infections - a critical step for malaria eradication.

Intrinsic multiplication rate variation and plasticity of human blood stage malaria parasites.

Pathogen multiplication rate is theoretically an important determinant of virulence, although often poorly understood and difficult to measure accurately. We show intrinsic asexual blood stage multiplication rate variation of the major human malaria parasite Plasmodium falciparum to be associated with blood-stage infection intensity in patients. A panel of clinical isolates from a highly endemic West African population was analysed repeatedly during five months of continuous laboratory culture, showing a range of exponential multiplication rates at all timepoints tested, mean rates increasing over time. All isolates had different genome sequences, many containing within-isolate diversity that decreased over time in culture, but increases in multiplication rates were not primarily attributable to genomic selection. New mutants, including premature stop codons emerging in a few isolates, did not attain sufficiently high frequencies to substantially affect overall multiplication rates. Significantly, multiplication rate variation among the isolates at each of the assayed culture timepoints robustly correlated with parasite levels seen in patients at clinical presentation, indicating innate parasite control of multiplication rate that contributes to virulence.

RecQ helicases in the malaria parasite Plasmodium falciparum affect genome stability, gene expression patterns and DNA replication dynamics.

The malaria parasite Plasmodium falciparum has evolved an unusual genome structure. The majority of the genome is relatively stable, with mutation rates similar to most eukaryotic species. However, some regions are very unstable with high recombination rates, driving the generation of new immune evasion-associated var genes. The molecular factors controlling the inconsistent stability of this genome are not known. Here we studied the roles of the two putative RecQ helicases in P. falciparum, PfBLM and PfWRN. When PfWRN was knocked down, recombination rates increased four-fold, generating chromosomal abnormalities, a high rate of chimeric var genes and many microindels, particularly in known 'fragile sites'. This is the first identification of a gene involved in suppressing recombination and maintaining genome stability in Plasmodium. By contrast, no change in mutation rate appeared when the second RecQ helicase, PfBLM, was mutated. At the transcriptional level, however, both helicases evidently modulate the transcription of large cohorts of genes, with several hundred genes-including a large proportion of vars-showing deregulated expression in each RecQ mutant. Aberrant processing of stalled replication forks is a possible mechanism underlying elevated mutation rates and this was assessed by measuring DNA replication dynamics in the RecQ mutant lines. Replication forks moved slowly and stalled at elevated rates in both mutants, confirming that RecQ helicases are required for efficient DNA replication. Overall, this work identifies the Plasmodium RecQ helicases as major players in DNA replication, antigenic diversification and genome stability in the most lethal human malaria parasite, with important implications for genome evolution in this pathogen.

Long read assemblies of geographically dispersed Plasmodium falciparum isolates reveal highly structured subtelomeres.

Background: Although thousands of clinical isolates of Plasmodium falciparum are being sequenced and analysed by short read technology, the data do not resolve the highly variable subtelomeric regions of the genomes that contain polymorphic gene families involved in immune evasion and pathogenesis. There is also no current standard definition of the boundaries of these variable subtelomeric regions. Methods: Using long-read sequence data (Pacific Biosciences SMRT technology), we assembled and annotated the genomes of 15 P. falciparum isolates, ten of which are newly cultured clinical isolates. We performed comparative analysis of the entire genome with particular emphasis on the subtelomeric regions and the internal var genes clusters.   Results: The nearly complete sequence of these 15 isolates has enabled us to define a highly conserved core genome, to delineate the boundaries of the subtelomeric regions, and to compare these across isolates. We found highly structured variable regions in the genome. Some exported gene families purportedly involved in release of merozoites show copy number variation. As an example of ongoing genome evolution, we found a novel CLAG gene in six isolates.  We also found a novel gene that was relatively enriched in the South East Asian isolates compared to those from Africa. Conclusions: These 15 manually curated new reference genome sequences with their nearly complete subtelomeric regions and fully assembled genes are an important new resource for the malaria research community. We report the overall conserved structure and pattern of important gene families and the more clearly defined subtelomeric regions.

Culture adaptation of malaria parasites selects for convergent loss-of-function mutants.

Cultured human pathogens may differ significantly from source populations. To investigate the genetic basis of laboratory adaptation in malaria parasites, clinical Plasmodium falciparum isolates were sampled from patients and cultured in vitro for up to three months. Genome sequence analysis was performed on multiple culture time point samples from six monoclonal isolates, and single nucleotide polymorphism (SNP) variants emerging over time were detected. Out of a total of five positively selected SNPs, four represented nonsense mutations resulting in stop codons, three of these in a single ApiAP2 transcription factor gene, and one in SRPK1. To survey further for nonsense mutants associated with culture, genome sequences of eleven long-term laboratory-adapted parasite strains were examined, revealing four independently acquired nonsense mutations in two other ApiAP2 genes, and five in Epac. No mutants of these genes exist in a large database of parasite sequences from uncultured clinical samples. This implicates putative master regulator genes in which multiple independent stop codon mutations have convergently led to culture adaptation, affecting most laboratory lines of P. falciparum. Understanding the adaptive processes should guide development of experimental models, which could include targeted gene disruption to adapt fastidious malaria parasite species to culture.

Extreme mutation bias and high AT content in Plasmodium falciparum.

For reasons that remain unknown, the Plasmodium falciparum genome has an exceptionally high AT content compared to other Plasmodium species and eukaryotes in general - nearly 80% in coding regions and approaching 90% in non-coding regions. Here, we examine how this phenomenon relates to genome-wide patterns of de novo mutation. Mutation accumulation experiments were performed by sequential cloning of six P. falciparum isolates growing in human erythrocytes in vitro for 4 years, with 279 clones sampled for whole genome sequencing at different time points. Genome sequence analysis of these samples revealed a significant excess of G:C to A:T transitions compared to other types of nucleotide substitution, which would naturally cause AT content to equilibrate close to the level seen across the P. falciparum reference genome (80.6% AT). These data also uncover an extremely high rate of small indel mutation relative to other species, primarily associated with repetitive AT-rich sequences, in addition to larger-scale structural rearrangements focused in antigen-coding var genes. In conclusion, high AT content in P. falciparum is driven by a systematic mutational bias and ultimately leads to an unusual level of microstructural plasticity, raising the question of whether this contributes to adaptive evolution.

Plasmodium falciparum clearance is rapid and pitting independent in immune Malian children treated with artesunate for malaria.

BACKGROUND: In Plasmodium falciparum-infected patients treated with artemisinins, parasitemia declines through so-called pitting, an innate splenic process that transforms infected red blood cells (iRBCs) into once-infected RBCs (O-iRBCs). METHODS: We measured pitting in 83 French travelers and 42 Malian children treated for malaria with artesunate. RESULTS: In travelers, O-iRBCs peaked at 107.7% initial parasitemia. In Malian children aged 1.5-4 years, O-iRBCs peaked at higher concentrations than in children aged 9-13 years (91.60% vs 31.95%; P = .0097). The parasite clearance time in older children was shorter than in younger children (P = .0001), and the decline in parasitemia in children aged 1.5-4 years often started 6 hours after treatment initiation, a lag phase generally absent in infants and older children. A 6-hour lag phase in artificial pitting of artesunate-exposed iRBCs was also observed in vitro. The proportion of iRBCs recognized by autologous immunoglobulin G (IgG) correlated with the parasite clearance time (r = -0.501; P = .0006) and peak O-iRBC concentration (r = -0.420; P = .0033). CONCLUSIONS: Antimalarial immunity correlates with fast artemisinin-induced parasite clearance and low pitting rates. In nonimmune populations, artemisinin-induced P. falciparum clearance is related to pitting and starts after a 6-hour lag phase. In immune populations, passively and naturally acquired immune mechanisms operating faster than pitting may exist. This mechanism may mitigate the emergence of artemisinin-resistant P. falciparum in Africa.

Generation of antigenic diversity in Plasmodium falciparum by structured rearrangement of Var genes during mitosis.

The most polymorphic gene family in P. falciparum is the ∼60 var genes distributed across parasite chromosomes, both in the subtelomeres and in internal regions. They encode hypervariable surface proteins known as P. falciparum erythrocyte membrane protein 1 (PfEMP1) that are critical for pathogenesis and immune evasion in Plasmodium falciparum. How var gene sequence diversity is generated is not currently completely understood. To address this, we constructed large clone trees and performed whole genome sequence analysis to study the generation of novel var gene sequences in asexually replicating parasites. While single nucleotide polymorphisms (SNPs) were scattered across the genome, structural variants (deletions, duplications, translocations) were focused in and around var genes, with considerable variation in frequency between strains. Analysis of more than 100 recombination events involving var exon 1 revealed that the average nucleotide sequence identity of two recombining exons was only 63% (range: 52.7-72.4%) yet the crossovers were error-free and occurred in such a way that the resulting sequence was in frame and domain architecture was preserved. Var exon 1, which encodes the immunologically exposed part of the protein, recombined in up to 0.2% of infected erythrocytes in vitro per life cycle. The high rate of var exon 1 recombination indicates that millions of new antigenic structures could potentially be generated each day in a single infected individual. We propose a model whereby var gene sequence polymorphism is mainly generated during the asexual part of the life cycle.

Optimized whole-genome amplification strategy for extremely AT-biased template.

Pathogen genome sequencing directly from clinical samples is quickly gaining importance in genetic and medical research studies. However, low DNA yield from blood-borne pathogens is often a limiting factor. The problem worsens in extremely base-biased genomes such as the AT-rich Plasmodium falciparum. We present a strategy for whole-genome amplification (WGA) of low-yield samples from P. falciparum prior to short-read sequencing. We have developed WGA conditions that incorporate tetramethylammonium chloride for improved amplification and coverage of AT-rich regions of the genome. We show that this method reduces amplification bias and chimera formation. Our data show that this method is suitable for as low as 10 pg input DNA, and offers the possibility of sequencing the parasite genome from small blood samples.

A subset of group A-like var genes encodes the malaria parasite ligands for binding to human brain endothelial cells.

Cerebral malaria is the most deadly manifestation of infection with Plasmodium falciparum. The pathology of cerebral malaria is characterized by the accumulation of infected erythrocytes (IEs) in the microvasculature of the brain caused by parasite adhesins on the surface of IEs binding to human receptors on microvascular endothelial cells. The parasite and host molecules involved in this interaction are unknown. We selected three P. falciparum strains (HB3, 3D7, and IT/FCR3) for binding to a human brain endothelial cell line (HBEC-5i). The whole transcriptome of isogenic pairs of selected and unselected parasites was analyzed using a variant surface antigen-supplemented microarray chip. After selection, the most highly and consistently up-regulated genes were a subset of group A-like var genes (HB3var3, 3D7_PFD0020c, ITvar7, and ITvar19) that showed 11- to >100-fold increased transcription levels. These var genes encode P. falciparum erythrocyte membrane protein (PfEMP)1 variants with distinct N-terminal domain types (domain cassette 8 or domain cassette 13). Antibodies to HB3var3 and PFD0020c recognized the surface of live IEs and blocked binding to HBEC-5i, thereby confirming the adhesive function of these variants. The clinical in vivo relevance of the HBEC-selected parasites was supported by significantly higher surface recognition of HBEC-selected parasites compared with unselected parasites by antibodies from young African children suffering cerebral malaria (Mann-Whitney test, P = 0.029) but not by antibodies from controls with uncomplicated malaria (Mann-Whitney test, P = 0.58). This work describes a binding phenotype for virulence-associated group A P. falciparum erythrocyte membrane protein 1 variants and identifies targets for interventions to treat or prevent cerebral malaria.

Induction of strain-transcending antibodies against Group A PfEMP1 surface antigens from virulent malaria parasites.

Sequence diversity in pathogen antigens is an obstacle to the development of interventions against many infectious diseases. In malaria caused by Plasmodium falciparum, the PfEMP1 family of variant surface antigens encoded by var genes are adhesion molecules that play a pivotal role in malaria pathogenesis and clinical disease. PfEMP1 is a major target of protective immunity, however, development of drugs or vaccines based on PfEMP1 is problematic due to extensive sequence diversity within the PfEMP1 family. Here we identified the PfEMP1 variants transcribed by P. falciparum strains selected for a virulence-associated adhesion phenotype (IgM-positive rosetting). The parasites transcribed a subset of Group A PfEMP1 variants characterised by an unusual PfEMP1 architecture and a distinct N-terminal domain (either DBLα1.5 or DBLα1.8 type). Antibodies raised in rabbits against the N-terminal domains showed functional activity (surface reactivity with live infected erythrocytes (IEs), rosette inhibition and induction of phagocytosis of IEs) down to low concentrations (<10 µg/ml of total IgG) against homologous parasites. Furthermore, the antibodies showed broad cross-reactivity against heterologous parasite strains with the same rosetting phenotype, including clinical isolates from four sub-Saharan African countries that showed surface reactivity with either DBLα1.5 antibodies (variant HB3var6) or DBLα1.8 antibodies (variant TM284var1). These data show that parasites with a virulence-associated adhesion phenotype share IE surface epitopes that can be targeted by strain-transcending antibodies to PfEMP1. The existence of shared surface epitopes amongst functionally similar disease-associated P. falciparum parasite isolates suggests that development of therapeutic interventions to prevent severe malaria is a realistic goal.

Selection of Plasmodium falciparum parasites for cytoadhesion to human brain endothelial cells.

Most human malaria deaths are caused by blood-stage Plasmodium falciparum parasites. Cerebral malaria, the most life-threatening complication of the disease, is characterised by an accumulation of Plasmodium falciparum infected red blood cells (iRBC) at pigmented trophozoite stage in the microvasculature of the brain(2-4). This microvessel obstruction (sequestration) leads to acidosis, hypoxia and harmful inflammatory cytokines (reviewed in (5)). Sequestration is also found in most microvascular tissues of the human body(2, 3). The mechanism by which iRBC attach to the blood vessel walls is still poorly understood. The immortalized Human Brain microvascular Endothelial Cell line (HBEC-5i) has been used as an in vitro model of the blood-brain barrier(6). However, Plasmodium falciparum iRBC attach only poorly to HBEC-5i in vitro, unlike the dense sequestration that occurs in cerebral malaria cases. We therefore developed a panning assay to select (enrich) various P. falciparum strains for adhesion to HBEC-5i in order to obtain populations of high-binding parasites, more representative of what occurs in vivo. A sample of a parasite culture (mixture of iRBC and uninfected RBC) at the pigmented trophozoite stage is washed and incubated on a layer of HBEC-5i grown on a Petri dish. After incubation, the dish is gently washed free from uRBC and unbound iRBC. Fresh uRBC are added to the few iRBC attached to HBEC-5i and incubated overnight. As schizont stage parasites burst, merozoites reinvade RBC and these ring stage parasites are harvested the following day. Parasites are cultured until enough material is obtained (typically 2 to 4 weeks) and a new round of selection can be performed. Depending on the P. falciparum strain, 4 to 7 rounds of selection are needed in order to get a population where most parasites bind to HBEC-5i. The binding phenotype is progressively lost after a few weeks, indicating a switch in variant surface antigen gene expression, thus regular selection on HBEC-5i is required to maintain the phenotype. In summary, we developed a selection assay rendering P. falciparum parasites a more "cerebral malaria adhesive" phenotype. We were able to select 3 out of 4 P. falciparum strains on HBEC-5i. This assay has also successfully been used to select parasites for binding to human dermal and pulmonary endothelial cells. Importantly, this method can be used to select tissue-specific parasite populations in order to identify candidate parasite ligands for binding to brain endothelium. Moreover, this assay can be used to screen for putative anti-sequestration drugs(7).

Design of a variant surface antigen-supplemented microarray chip for whole transcriptome analysis of multiple Plasmodium falciparum cytoadherent strains, and identification of strain-transcendent rif and stevor genes.

BACKGROUND: The cytoadherence of Plasmodium falciparum is thought to be mediated by variant surface antigens (VSA), encoded by var, rif, stevor and pfmc-2tm genes. The last three families have rarely been studied in the context of cytoadherence. As most VSA genes are unique, the variability among sequences has impeded the functional study of VSA across different P. falciparum strains. However, many P. falciparum genomes have recently been sequenced, allowing the development of specific microarray probes for each VSA gene. METHODS: All VSA sequences from the HB3, Dd2 and IT/FCR3 genomes were extracted using HMMer software. Oligonucleotide probes were designed with OligoRankPick and added to the 3D7-based microarray chip. As a proof of concept, IT/R29 parasites were selected for and against rosette formation and the transcriptomes of isogenic rosetting and non-rosetting parasites were compared by microarray. RESULTS: From each parasite strain 50-56 var genes, 125-132 rif genes, 26-33 stevor genes and 3-8 pfmc-2tm genes were identified. Bioinformatic analysis of the new VSA sequences showed that 13 rif genes and five stevor genes were well-conserved across at least three strains (83-100% amino acid identity). The ability of the VSA-supplemented microarray chip to detect cytoadherence-related genes was assessed using P. falciparum clone IT/R29, in which rosetting is known to be mediated by PfEMP1 encoded by ITvar9. Whole transcriptome analysis showed that the most highly up-regulated gene in rosetting parasites was ITvar9 (19 to 429-fold up-regulated over six time points). Only one rif gene (IT4rifA_042) was up-regulated by more than four fold (five fold at 12 hours post-invasion), and no stevor or pfmc-2tm genes were up-regulated by more than two fold. 377 non-VSA genes were differentially expressed by three fold or more in rosetting parasites, although none was as markedly or consistently up-regulated as ITvar9. CONCLUSIONS: Probes for the VSA of newly sequenced P. falciparum strains can be added to the 3D7-based microarray chip, allowing the analysis of the entire transcriptome of multiple strains. For the rosetting clone IT/R29, the striking transcriptional upregulation of ITvar9 was confirmed, and the data did not support the involvement of other VSA families in rosette formation.